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Abstract 

The time evolution of cosmological parameters in early Universe at the deconfinement transition is studied 
by an equation of state (EoS) which takes into account the finite baryon density and the background magnetic 
field. The non perturbative dynamics is described by the Field Correlator Method (FCM) which gives, with 
a small number of free parameters, a good fit of lattice data. The entire system has two components, i.e. the 
quark-gluon plasma and the electroweak sector, and the solutions of the Friedmann equation show that the 
scale factor, a(t), and H{t) = {l/a)da/dt are weakly dependent on the EoS, but the deceleration parameter, 
q{t), and the jerk, j{t), are strongly modified above the critical temperature T)., corresponding to a critical 
time tc — 20 — 25p,s. The time evolution of the cosmological parameters suggest that above and around Tc 
there is a transient state of acceleration typical of a matter dominated Universe; this is entailed by the QCD 
strong interaction driven by the presence of massive colored objects. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Nq 


INTRODUCTION 

Lattice simulations of Quantum Chromody¬ 
namics (QCD) and phenomenological analyses of 
relativistic heavy ion collisions data clearly indi¬ 
cate that the transition from quarks and gluons 
to colorless states occurs in a non perturbative 
regime such that for large baryon chemical po¬ 
tential, p,B, the transition is first order and for 
small pLB it is cross-over at a (pseudo) critical 
temperature U ~ 155 ± 15 Mev [H, 1^. 

Moreover lattice QCD in a background mag¬ 
netic field , B, shows [ 1 , that the transition 
temperature is reduced by the magnetic held and 
that the thermal QCD medium is paramagnetic 
around and above the transition temperature. 

The previous parameters, T,piB and B, give 
different and important informations on the dy¬ 
namics of the phase transition which, in turn, 
has strong implications at cosmological level in 
the early Universe ( see j^,|^ for a recent review). 

In particular, the equation of state (EoS) in 
the quark-gluon plasma phase affects directly 
the density fluctuactions of the thermodynamic 
quantities 0-0 and the emission of gravita¬ 


tional waves 0 , [ll- 15] during the cosmological 
evolution. 

In ref. 0 the density huctuactions in the 
early Universe have been studied by using the 
MIT bag model EoS and the old lattice QCD 
EoS with a first order deconfinement phase tran¬ 
sition. A more recent and realistic QCD EoS 
has been introduced in ref. 0 improving the re¬ 
liability of the calculations of the fluctuactions 
and of the modihcations in the time evolution of 
the cosmological scale factor a{t). The effects of 
the external magnetic field on the energy density 
fluctuations have been considered in ref. ibya 
phenomenological QCD EoS which includes B. 
An analogous analysis for finite density ^b with 
different variants of MIT bag model has been 
carried out in ref. 


m- 


The results of the previous studies show a 
smooth time dependence of the thermodynamic 
quantities, however they essentially consider the 
effect of the transition on the thermodynamic 
fluctuactions and on the scale factor. 


In this paper we shall analize the modihca¬ 
tions on the early Universe evolution, by using 
a realistic EoS which depends on T,pLB and B, 
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with a particular focus on the evolution of the 
cosmological parameters rather than on the fluc- 
tuactions. 

To describe the EoS in the quark-gluon 
plasma phase we shall apply the Field Corre¬ 
lator Method (FCM) [I^. This specific choice is 
motivated by the non perturbative dynamics of 
the FCM which describes, with a small number 
of free parameters, lattice data of the thermody¬ 
namics quantities ( pressure and energy density) 
at finite temperature and their dependence on 


1. COSMOLOGICAL PARAMETERS 

The parameters H{t) = a /a, g'(t)and j{t) 
can naturally be defined making use of Taylor 
series for the scale factor a{t) near a generic time 
t*: 

a{t) = a{t*) + a — t*) + 

+ + ■•■ ( 3 ) 


fiB and B [13-Il9l]. Here we willalso show that ^hich can be written as 


its extension to a finite B 20|, l2]| capture cor¬ 


rectly the main features of IQCD EoS under a 
magnetic field. 

The discussion of the combined effects of the 
chemical potential and of the external magnetic 
field, during the deconfinement transition, in the 
Friedmann equations is of interest in its own 
right. However, as we shall see, the most in¬ 
teresting results are on the behavior of the de¬ 
celeration parameters, q, and the jerk, j, defined 
as {' indicates the time derivative) 


q = 


a a 

n ^ 


a a 


J = 


( 1 ) 


( 2 ) 


which have an important role in describing 
the cosmological evolution 
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Indeed, it turns out that the deceleration and 
the jerk are strongly modified above the critical 
temperature Tc, corresponding to a critical time 
tc — 20 — 25/is, and that the EoS and the time 
evolution of the cosmological parameters suggest 
that above Tc the system has the typical behav¬ 
ior of a matter dominated Universe, with clusters 
of colored particles. 

The plane of the paper is as follows. In Sec. 1 
we shall consider the relevant cosmological equa¬ 
tions and parameters. The FCM, with the de¬ 
pendence on fiB and B is recalled and compared 
with lattice data in Sec.2. Sec.3 is devoted to 
the solution of the Friedmann equation with the 
FCM EoS and Sec.4 contains the electroweak 
contributions to the EoS of the entire dynamical 
system. In Sec.5 the results on the cosmological 
parameters are discussed and our comments and 
conclusions are in Sec.6. 


+ 4{ji?’)(i-)((-Cf+ ... (4) 

Basic characteristics of the cosmological evolu¬ 
tion, both static and dynamical, can be ex¬ 
pressed in terms of Hq, the present time value 
of H{t), and of the deceleration qq. The other 
parameters, i.e. the higher time derivatives 
of the scale factor, enable to construct model- 
independent kinematics of the cosmological ex¬ 
pansion. 

Indeed cosmological models can be tested 
by expressing the Friedmann equation in terms 
of directly measurable cosmological scalars con¬ 
structed out of higher derivatives of the scale 
factor, i.e q,j,.. 


mM 


To illustrate this aspect, let us consider a 
simple two-component Universe filled with non- 
relativistic matter with density M^/a^ and ra¬ 
diation with density Mr ja^ ( Mm , Mr constants) 
which do not interact with each other 2^ and 
without a magnetic background. By writing the 
Friedmann equation in the form (SvrG/S = 1) 


a ^ k Mm Mr 
a? ^ a? a'^ 


( 5 ) 


and diffrentiating twice with respect to t, one 
gets 


1 Mm Mr 


a =- 


and 


2 


m Mjyi f Mr / 

CL = -“h 3— -^d . 


( 6 ) 


( 7 ) 
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Then the deceleration and the jerk can be writ¬ 
ten as q = Ajl + R and j = A + 3R where 
A = and R = Mr/a^H"^. For a flat 

Universe filled only with non relativistic matter 
one has R = 0, q = 1/2, A = 1 and j = 1; if 
one considers only radiation,then ^ = 0, g = 1, 
R = 1 and j = 3. A deviation from these val¬ 
ues of q and j indicates a mixture of matter and 
radiation and/or an interaction between the two 
components.Indeed the same values can be eas¬ 
ily obtained by a simple application of the Fried¬ 
mann equation (for k = 0), written as 


a' de/dt /SvrG ^ 

7 ^"3(e + p) ^ 


( 8 ) 


by the EoS p = c^e with a constant speed of 
sound. The analitic solution is 


e{t) 


_4e(to)_ 

[3y^^^(l + c2)(t-to) + 2]' 


( 9 ) 


and 


a{t) r3 
~ 1-9 


87rGe(to) 


a(to) '■2 
and one obtains 

q = 


{l + cl){t-to) + l] 3(1+'=*) 
( 10 ) 


o"a 




l + 3ci 


( 11 ) 


a"'a2 


J = 


^/3 


= ( k ^){2 + 3 a ( 12 ) 


which reproduce the previous values for matter 
dominated (c^ = 0) and radiation dominated 
(Cg = 1/3) Universe. 

In general the speed of sound is not constant 
and in the next sections we shall discuss the be¬ 
havior of the cosmological parameters q and j 
during the deconfinement transition on the ba¬ 
sis of the Friedmann equation (8) and by using 
the energy density, e and the pressure, p, in the 
FCM, after fitting the QCD lattice data at finite 
temperature, pB and B. 


2.riELD CORRELATOR METHOD 

Many phenomenological models of QCD at 
finite temperature and density cannot make re¬ 
liable predictions for the two relevant limits, i.e. 


high temperature and small chemical potential 
or high chemical potential and low tempera¬ 
ture. This is clearly a serious drawback, since 
those models cannot be fully tested. One of the 
few exceptions is the Field Correlator Method 
(FCM) [3], which is able to cover the full 
temperature-chemical potential-magnetic back¬ 
ground field space and contains ab initio the 
property of confinement, which is expected to 
play a role, at variance with other models like, 
e.g., the Nambu - Jona Lasinio model. 

Indeed, the approach based on the FCM pro¬ 
vides a systematic tool to treat non perturba¬ 
tive effects in QCD by gauge invariant field cor¬ 
relators and gives a natural treatment of the 
dynamics of confinement (and of the deconfine¬ 
ment transition) in terms of the Gaussian, i.e. 
quadratic in the tensor correlators for the 
chromo-electric (CE) field, and Df, and for 
the chromo-magnetic field (CM), and . 

In particular, these correlators are related to 
the simplest non trivial 2-point correlators for 
the CE and CM fields by 


< Trf[Ci{x)^{x,y)Ck{y)^{y,x)] >= 


{z) + Di (z) + Z'l 


C/ 3 I 


± ZiZk 


dz^ 

dD^{z) 

dz^ 


(13) 


where z = x — y and C indicates the CE (E) 
field or CM (H) field (the minus sign in the previ¬ 
ous expression corresponds to the magnetic case) 
and 


ry 

= Pexp[ig / A'^dZfj,] (14) 


is the parallel transporter. 

The FCM has been extended to finite temper¬ 
ature and chemical potential 17Hl9l| and the an¬ 
alytical results, in the gaussian approximation, 
are in good agreement with the lattice data on 
thermodynamic quantities, (available for small 
Pb only). Moreover, the application of the FCM 
for large values of the chemical potential allows 
to obtain a simple expression of the Equation of 
State of the quark- gluon matter in the range 
of the baryon de nsity relevant for the study of 

[iF 


neutron stars 


2a, 
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pressure, p, evaluated by the FCM. speed of sound evaluated by the FCM. 



FIG. 2: Comparison with lattice data [lilii of 
the interaction measure A = (e — 3p)/T* evaluated 
by the FCM. 


The comparisons with lattice data of the 
FCM predictions for the pressure (p/r^),the in¬ 
teraction measure, A/T^ = (e —3p)/T^, and the 
speed of sound, at = 0 and /x_b = 0.4GeV, 
are depicted in figs. 1,2 and 3. 

More recently, the effect of a background 


magnetic field, B, has been included in the FCM 
equation of state 0, 0] and the quark (q) pres¬ 
sure and the gluon (g) pressure turn out to be 


Pqi^) = 

2 (A(/x) -b A(-/x)) 

3 eqB 

+'r(-(j))]. (15) 




I 


z^dz 


3-^ Jo exp(. + |^)- 


1 


(16) 

where </>(/x) ,A(/x) and r(/i) are respectively given 
by 

Pzdpz 


= [ 

Jo 

Jo 


0 l-bexp(£^) 

1 


(17) 


0 . /p2 _|_ fji2 


q exp 




+ 1 
(18) 


PC 

t{p) = 

Jo 


dpz 




+ 1 
(19) 
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FIG. 4: Comparison of the FCM calculations for the 
entropy density with lattice data 0 for QCD in a 
background magnetic field B 


with /i = /i — , 


= TTig + egB with rUn 


the 


vviui.1 utq 

current mass {nid = 5GeV, ruu = lOMeV and 
ms = 140MeV), and 


Vi{T) = co + f{T/Tg) 

= 0.15 + 0.175 ^1.35^-!^ (20) 

is the quark-quark interaction potential. 

In the limit eB —)• 0 the only non zero terms 
are A(/r) and A(—/x) and t he p revious formulas 
reproduce the case ur 0 [l7|. 

In fig. 4 one compares the ratio s/T^ (s be¬ 
ing the entropy density) evaluated in the FCM 
with lattice data for different values of the back¬ 
ground field (iil . 

We consider a value of eB = O.GGeV^ that 
corresponds to a maximum estimate of the mag¬ 
netic field generated at the early times [^. In 
some possible scenario also a finite could be 
still significant ~ T) at the QCD transition 
@]. We consider a = 0.4 GeV, but we will 
see that its impact is however quite limited. 

The deconfinement temperature depends on 
/X B and B and therefore the critical time tg of the 
transition turns out to be 24 ^s, 21.5 /xs, 15.6 /xs 
and 13.6/xs respectively for ^b = O.OGeV and 


B = 0.0, /xb = 0.4GeV and B = 0.0, /xb = 0.0 
and B = O.GGeV^, and /xb = 0.4 GeV and B = 
0.6 GeV^. 


3.FRIEDMANN EQUATION AND QUARK- 
GLUON PLASMA IN THE FCM 

In the FCM the EoS of the quark-gluon 
plasma, i.e. p(e), is obtained by direct calcu¬ 
lations of p{T) and e{T), which inserted in the 
Friedman equation (8) give the time dependence 
of the temperature , T{t), and the corresponding 
time evolution of the thermodynamic quantities. 

For the initial conditions ti = 1/xs, Ti ~ 500 
MeV and e* ~ 110 GeV/fm^, the function T{t), 
solution of the Friedman equation, is shown in 
fig.5 for different values of /xb and compared 
with the MIT bag model with the same ini¬ 
tial conditions and a bag pressure Bmu = 220 
MeV/fmV 

The arrows in fig.5 (and in the next fig¬ 
ures) correspond to the critical temperatures 
for the different specific sets of the parame¬ 
ters and therefore to the corresponding criti¬ 
cal time, tg, when the phase transition occurs. 
The critical temperature Tg = 160 MeV, for 
/Xb = 0,corresponds to tg ~ 25/xs, which de¬ 
creases to tg ~ 22/xs for /XB = 0.4 GeV. The 
curves are also plotted for t > tg, i.e. for tem¬ 
perature below the transition point although the 
effective degrees of freedom below Tg are non in¬ 
cluded in the present paper. 

The equation of state p/e in the FCM, re¬ 
ported in fig. 6, shows a small dependence on 
/Xb up to 400 MeV. 

The role of the background magnetic field in 
the calculation of T{t) and in the time evolu¬ 
tion of the EoS are reported in fig.7 and in fig.8, 
which also describe the combined effect of the 
finite density and of the background magnetic 
field. 

Notice that the results in the MIT bag model 
( green curves in figs. 6,7,8) are clearly different 
from the EoS evaluated by the ECM and that the 
magnetic field B can move down tg by nearly a 
factor of 2. 
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FIG. 5: T{t) solution of the Friedmannn equation (8) 
by using the EoS in the FCM. The arrows indicate 
the critical time, tc of the transition given in the text. 



FIG. 6: EoS in the EGM for finite chemical poten¬ 
tial. The green curve referes to the MIT bag model 
calculation with the parameters given in the text. 


4.ELECTROWEAK CONTRIBUTION TO 
THE EQUATION OF STATE 

The time evolution of the cosmological pa¬ 
rameters depend on the EoS of the entire sys- 



FIG. 7: The temperature profile T{t) for different 
values of /is and B. 



FIG. 8: Time evolution of the EoS in the FCM for 
finite chemical potential and background magnetic 
field. 

terns and therefore one has to take into account 
not only the quark-gluon plasma (qgp) degrees 
of freedom, discussed in the previous section, but 
also the contributions of the electroweak sector 
to the pressure, Pew and the energy density, €ew 
The total pressure and energy density are there- 
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FIG. 9: T{t) with/without the electroweak contribu- FIG. 10: EoS with/without the electroweak contribu¬ 
tion and the background magnetic field. tion and the background magnetic field, for = 0 


fore ptot — Pggp Pew (^tot — ^ggp ^ew’ 

The electroweak sector is described as a sys¬ 
tem of free massless particles, i.e. 


^ew 


TT 


Qew - T 

yew 


( 21 ) 


Pew — 9ew , (22) 

with the number of degrees of freedom pew = 
14.45. 

The temperature prohle is almost unmodified 
by the introduction of the electroweak contribu¬ 
tions as one can see from fig.9, where the curves 
are plotted for = 0 , since the chemical poten¬ 
tial produces negligible changes (see figs. 5,6). 

On the other hand, the EoS of the entire sys¬ 
tem strongly reflects the presence of the elec¬ 
troweak terms. In fig. 10 the time evolution of 
p/e for the entire system is depicted and, indeed, 
one immediately notes that in the time interval 
t ~ 10 — 25ps the electroweak part has a minor 
role, but near the critical temperature there is a 
clear change in the shape and the EoS tends to 
p = e/3, i.e. to a radiation dominated Universe, 
for long time. The curves are plotted for time 
longer than tc, just to show the behavior below 


Tc without the contribution of the colorless effec¬ 
tive degrees of freedom after the transition, not 
included in the present analysis and that will be 
discussed in a forthcoming paper. 


5.EVOLUTION OF THE COSMOLOGICAL 
PARAMETERS 

The evolution of the cosmological parameters 
during the deconfinement transition is directly 
related with the EoS. In fig. 11 and in fig. 12 
are respectively depicted the time behavior of 
the scale factor and of H(t) for different values of 
the chemical potential, of the magnetic field and 
also with/without the electroweak contribution. 
The hnal result is essentially independent on the 
specific setting. 

However the deceleration q{t) and the jerk 
j{t) strongly follows the time evolution of the 
EoS (in hg. 10). Both q and j ( see fig. 13 
and hg. 14) after initial values corresponding to 
a radiation dominated Universe ((7 ~ 1, j ~ 3), 
tend to q ~ 1/2 and j ~ 1 , typical of a mat¬ 
ter dominated Universe, approaching from above 
the critical temperature. Indeed, without the 
electroweak sector the values of the cosmological 
parameters for time t ~ 100/is would be the typ- 
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FIG. 11: Time evolution of the scale factor for finite 
chemical potential, with/without the electroweak 
contribution and the background magnetic field. 


FIG. 13: q{t) for finite chemical potential, 

with/without the electroweak contribution and the 
background magnetic field. 



FIG. 12: H{t) for finite chemical potential, 
with/without the electroweak contribution and the 
background magnetic field. 



FIG. 14: j{t) for finite chemical potential, 

with/without the electroweak contribution and the 
background magnetic field. 


ical ones of a matter dominated Universe. How¬ 
ever , near the transition point the electroweak 
terms in the EoS start to be relevant and there¬ 
fore q ^ 1 and J —?• 3 for long time. 

A possible interpretation of this peculiar be¬ 


havior can be done with the help of the sim¬ 
ple model in Sec.l, where A = —2j + 6q and 
R = j — 2q are respectively the matter and ra¬ 
diation terms in the Friedmann equation. For a 
matter dominated flat Universe A = 1 and R = 0 
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FIG. 15: Time evolution of A and H according to the 
simple model in eq.(5) but using q(t) and j(t) in the 
FCM. 


whereas for a radiation dominated one A = 0 
and R = 1 and since in eq.(5) there is no inter¬ 
action between the two components A + R = 1. 
Of course, this condition in not satisfied by us¬ 
ing q and j calculated in the FCM and different 
values of A and R signal a mixture and/or an 
interaction between the two components. 

By using the values of q and j in the FCM, 
the time evolution of A and R is given in fig. 15. 
R decreases from an initial value ~ 1 to ~ 0.3 
at t ~ 25/US and to ~ 0 at t ~ 32/is and then 
quickly reaches again the value of a radiation 
dominated Universe. A has the corresponding 
evolution , starting from A ~ 0.0, increasing to a 
value ~ 1 at t ~ 25/us, with a maximum ^ ~ 1.8 
at t ~ 35/US, and finally decreasing to ~ 0,i.e. 
a radiation dominated Universe. 

Therefore there is a clear mixture of the two 
components which strongly interact each other. 

The observed behavior of A and R implies 
that above the transition and before the domi¬ 
nance of the electroweak sector, the system has 
essentially the EoS of interacting matter. On 
the other hand, above Tc the color degrees of 
freedom are still not neutralized and therefore 
the previous results suggest the formation of col¬ 
ored and massive clusters near the deconfine¬ 


ment transition before the formation of colorless 
bound states. 

Indeed, this is a well known interpretation 
of the QCD EoS at finite temperature in terms 
of quasi-particles where quasi-quarks and quasi¬ 
gluons have a dynamical, temperature depen¬ 
dent, effective masses which mimic the interac¬ 
tion and that near Tc are large, i.e. in the range 
~ 0.6- 1.2 GeV for = 0 = 5 [HI, [H, 0 [35 1. 


COMMENTS AND CONCLUSIONS 


The analysis in the previous sections shows 
that during the deconfinement transition the 
time evolution of the scale factor and of H (t) are 
weakly sensitive to the EoS and that , on the 
contrary, the cosmological parameters q{t) and 
j{t) follows the behavior of the ratio p/e. Start¬ 
ing from a radiation dominated Universe, the 
time evolution of the EoS indicates that above 
and near the transition time tc the entire system 
(quark gluon plasma -|- electroweak sector) is in a 
matter dominated state (g « 0.5 and j « 1). For 
longer time, the evolution is again dominated by 
the radiation EoS. 

The introduction of a finite baryon chemical 
potential and a background magnetic field do not 
qualitatively change this dynamical picture. 

On the other hand, since above and near 
Tc one has color degrees of freedom, the mat¬ 
ter state which drive the EoS is , presumably, 
formed by color massive objects, as suggested 


by the quasi-particle models |3ll-l33l|. Indeed, 


the behavior of A{t), i.e. of the matter content 
of the system, reported in fig. 15 is analogous 
to the time evolution of the interaction measure 
(e — 3p)/T^, in fig. 16, where the electroweak 
sector has no role. 

The effective degrees of freedom below Tc, 
i.e. for t > tc, are not included in the present 
study and the extrapolations for long time give 
the time evolution of a radiation dominated Uni¬ 
verse. On the other hand, below Tc and at small 
baryon chemical potential one expects the for¬ 
mation of qq states which decay in electroweak 
final states. It is known that a hadron resonance 
gas model is able to give a good desertion of the 
QCD matter down to T ~ l/3rc 0, 37|. Includ- 
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FIG. 16: Time evolution of the interaction measure. 

ing properly the EoS of hadronic matter proba¬ 
bly will lead to further extend the stage of matter 
dominated dynamics even beyond t ~ lOOfm/c. 
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